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We study the fractal and multifractal properties (i.e. the generalized dimensions of the harmonic 
measure) of a 2-parameter family of growth patterns that result from a growth model that interpo- 
lates between Diffusion Limited Aggregation (DLA) and Laplacian Growth Patterns in 2-dimensions. 
The two parameters are /3 which determines the size of particles accreted to the interface, and C 
which measures the degree of coverage of the interface by each layer accreted to the growth pattern 
at every growth step. DLA and Laplacian Growth are obtained at /3 = 0, C = and P = 2, C — 1, 
respectively. The main purpose of this paper is to show that there exists a line in the /3 — C phase di- 
agram that separates fractal (D < 2) from non-fractal (D=2) growth patterns. Moreover, Laplacian 
Growth is argued to lie in the non-fractal part of the phase diagram. Some of our arguments are not 
rigorous, but together with the numerics they indicate this result rather strongly. We first consider 
the family of models obtained for /3 = 0, C > 0, and derive for them a scaling relation D — 2Dz. 
We then propose that this family has growth patterns for which D = 2 for some C > Ca, where 
Ccr may be zero. Next we consider the whole (3 — C phase diagram and define a line that separates 
2-dimensional growth patterns from fractal patterns with D < 2. We explain that Laplacian Growth 
lies in the region belonging to 2-dimensional growth patterns, motivating the main conjecture of 
this paper, i.e. that Laplacian Growth patterns are 2-dimensional. The meaning of this result is 
that the branches of Laplacian Growth patterns have finite (and growing) area on scales much larger 
than any ultra-violet cut-off length. 



I. INTRODUCTION 

In recent work 1^, ^ we have introduced a model of 
fractal growth processes that interpolates between Dif- 
fusion Limited Aggregation (DLA) [|j and Laplacian 
Growth Patterns || , and employed this model to show 
that these processes are not in the same universality 
classes. The aim of this paper is to study the fractal 
properties of the resulting clusters. In particular we will 
be led to conjecture that Laplacian Growth is asymptot- 
ically of dimension 2, and in this sense is not a fractal 
at all. This is in contradistinction to DLA for which the 
dimension had been computed to be 1.713...[|j. 

Laplacian Growth Patterns are obtained when the 
boundary F of a 2-dimensional domain is grown at a 
rate proportional to the gradient of a Laplacian field P. 
Outside the domain V^P = 0, and each point of P is 
advanced at a rate proportional to VP ||, ||]- In Dif- 
fusion Limited Aggregation (DLA) Q a 2-dimensional 
cluster is grown by releasing fixed size random walkers 
from infinity, allowing them to walk around until they 
hit any particle belonging to the cluster. Since the par- 
ticles are released one by one and may take arbitrarily 
long time to hit the cluster, the probability field is quasi- 
stationary and in the complement of the cluster we have 
again V^P = 0. The boundary condition at infinity is the 
same for the two problems; in radial geometry as r — > oo 
the flux is VP = const x r/r. Since the probability for 
a random walker to hit the boundary is again propor- 
tional to |VP|, one could think that in the asymptotic 
limit when the size of the particle is much smaller than 
the radius of the cluster, repeated growth events lead to 
a growth process which is similar to Laplacian Growth. 



Of course, the ultraviolet regularizations in the two pro- 
cesses were taken different; in studying Laplacian Growth 
one usually solves the problem with the boundary con- 
dition P = an where a is the surface tension and k the 
local curvature of P Q . Without this (or some other) ul- 
traviolet regularization Laplacian Growth reaches a sin- 
gularity (cusps) in finite time Q . In DLA the ultraviolet 
regularization is provided by the finite size of the random 
walkers. However, many researchers believed Q that this 
difference, which for very large clusters controls only the 
smallest scales of the fractal patterns, were not relevant, 
expecting the two models to lead to the clusters with the 
same asymptotic dimensions. While we argued recently 
that the difference in ultraviolet regularization is indeed 
not crucial the two problems are nevertheless in two 
different universality classes. To establish this we have 
constructed a family of growth processes that includes 
DLA and a discrete version of Laplacian Growth as ex- 
treme members, using the same ultraviolet regularization 
(and see Sect. 2 for a further discussion of the regulariza- 
tion). We thus exposed the essential difference between 
DLA and Laplacian Growth. DLA is grown serially, with 
the field being updated after each particle growth. On 
the other hand all boundary points of a Laplacian pat- 
tern are advanced in parallel at once (proportional to 
VP). We showed that this difference is fundamental to 
the asymptotic dimension, putting the two problems in 
different universality classes [Q . Here we wish to go fur- 
ther and suggest that Laplacian Growth patterns are 2- 
dimensional. 

In Sect. 2 we review briefly the two parameter model 
that had been introduced to establish these results. We 
discuss there the two parameters (3 and C that are used 
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to interpolate between DLA and Laplacian Growth. In 
Sect. 3 we analyze the generalized dimensions Dg and re- 
late them to the scaling of moments of objects which are 
natural to the theory. In Sect. 4 we discuss first a family 
of growth models which is a 1-parameter generalization 
of DLA, (/3 = 0, < C < 1, and show that the fractahty 
of DLA is lost for some C > Ca in favor of 2-dimensional 
growth patterns. It is not impossible that Ccr = 0. For 
growth patterns in this family we derive a scaling relation 
D = 2D3. Under some plausible assumptions we propose 
that for C > Ccr there exists another scaling relation, i.e. 
D = 1 + D2, which implies immediatly that D = 2. Sec- 
ondly we discuss the 1-parameter family of models that 
generalizes Laplacian Growth {(3 — 2, < C < 1 and 
show that the above relation is not obtained here, lead- 
ing to the existence of fractal patterns also for high values 
of C. Finally in Sect. 5 we reach the main conjecture of 
this paper, i.e. that Laplacian Growth patterns are 2- 
dimensional. In Sect. 6 we offer a discussion and some 
open questions that are left for future research. 



II. ITERATED CONFORMAL MAPS FOR 
PARALLEL GROWTH PROCESSES 

The method of iterated conformal maps for DLA was 
introduced in j|] . In j^, ^ we have presented a generaliza- 
tion to parallel growth processes. We were interested in 
(w) which conformally maps the exterior of the unit 
circle e** in the mathematical w-plane onto the comple- 
ment of the (simply-connected) cluster of n particles in 
the physical z-plane. The unit circle is mapped onto the 
boundary of the cluster. In what follows we use the fact 
that the gradient of the Laplacian field VP(z(s)) is 



|VP(z(s))| 



|$(n)'(e«e)| 



z(s) = $(")(£*'') 



(1) 



Here s is an arc-length parametrization of the boundary. 
The map $(")(z«) is constructed recursively. Suppose 
that we have already <I>'"^(ti;) which maps to the exterior 
of a cluster of n particles in the physical plane and we 
want to find the map <I>("+p) (ui) after p additional par- 
ticles were added to its boundary at once, each propor- 
tional in size to the local value of \\/Pf/^. To grow one 
such particle we employ the elementary map (px.e which 
transforms the unit circle to a circle with a semi-spherical 
"bump" of linear size \/A around the point w = e^^: 
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(2) 
(3) 



If we update the field after the addition of this single 
particle, then 



$("+l)(^,)=ci>(»)(0;,„^^,,„^^(u;)) 



(4) 



where $(")(e*^"+i) is the point on which the (n -I- l)-th 
particle is grown and ^/X^ is the size of the grown particle 
divided by the Jacobian of the map, $ (e'^"+i), at that 
point. 

The map $'^"+^-'(w) adds on a new semi-circular bump 
to the image of the unit circle under (w). The 
bumps in the z-plane simulate the accreted particles 
in the physical space formulation of the growth pro- 
cess. For the height of the bump to be proportional to 
I VP(z(s))|''/2 we need to choose its area proportional to 
|$(«-i)'(e<e„)|-/3 (seeEq. (0)), or 



Xn 



An 



|$(n-l)'(gie„)|/3+2 



(5) 



With P — these rules produce a DLA cluster, for which 
the particles are of constant area. With /9 = 2 we grow 
bumps in the physical space whose linear scale is pro- 
portional the gradient of the field \W P{z{s))\, as is ap- 
propriate for Laplacian Growth. Next, to grow p (non- 
overlapping) particles in parallel, we accrete them with- 
out updating the conformal map. In other words, to add 
a new layer of p particles when the cluster contains m 
particles, we need to choose p angles on the unit circle 
{(^m+k}k=i- At these angles we grow bumps which in the 
physical space have the wanted linear scale (ranging from 
constant to proportional to the gradient of the field): 



A 



An 



|$(m) (e'''" + fc)|/3+2 



1,2. 



(6) 



After the p particles were added, the conformal map and 
thus the field should be updated. In updating, we will 
use p compositions of the elementary map 4)\fi{w). Of 
course, every composition effects a reparametrization of 
the unit circle, which has to be taken into account. To 
do this, we define a series {9m+kYk^i according to 



(7) 



Next we define the conformal map used in the next layer 
growth according to 

<i>(™+P)(c^) S $(™)O0,„^,.;,„^,O. . .0</>,,„^^,,„^^(c^) . (8) 

In this way we achieve the growth at the images under 
$(™) of the points {dm+k\^k=i- compute the 9 series 
from a given 9 series we use Eq.(|^) to rewrite Eq.(|^) in 
the form 



„A„+,(e^''"+'=) (9) 



The inverse map 0g \ is given by 0g ^(w) = 
e'Voi(e"'*w) with 



Aa;2 ± ^A2w4 - cj2[i - (l + A)cj2][a;2 - (1 + A)] 



1-(1 + A)t^2 



(10) 
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where the positive root is taken for Re w > and the 
negative root for Re w < 0. 

Evidently, Laplacian Growth calls for choosing the se- 
ries {drn+k}k=i such as to havc full coverage of the unit 
circle (implying the same for the boundary T). On the 
other hand DLA calls for growing a single particle before 
updating the field. Since it was shown that in DLA 
growth A„ decreases on the average when n increases, in 
the limit of large clusters DLA is consistent with van- 
ishingly small coverage of the unit circle. To interpolate 
between these two cases we introduce a parameter that 
serves to distinguish one growth model from the other, 
giving us a 2-parameter control (the other parameter is 
This parameter is the degree of coverage. Since the 
area covered by the pre-image of the n-th particle on 
the unit circle is approximately 2y/X^, we introduce the 
parameter 



m-\-k 



(11) 



fe=i 



(In we showed how to measure the coverage exactly). 
Since this is the fraction of the unit circle which is covered 
in each layer, the limit of Laplacian Growth is obtained 
with C — 1. DLA is asymptotically consistent with C — 0. 
Of course, the two models differ also in the size of the 
growing bumps, with DLA having fixed size particles, 
(/3 = in Eq.(|)), and Laplacian Growth having particles 
proportional to VP (/3 = 2 in Eq.(||)). Together with C 
we have a two parameter control on the parallel growth 
dynamics, with DLA and Laplacian Growth occupying 
two corners of the (3,C plane, at the points (0,0) and 
(2,1) respectively. 

Obviously, the partially serial growth within the layer 
introduces an additional freedom which is the order of 
placement of the bumps on the unit circle. In jl], |^ we 
have shown that the order is in fact immaterial as far 
as the asymptotic fractal properties of the clusters are 
concerned. Accordingly, we will take random choices of 
(^m+k with a rule of skipping overlaps. 

We should note that in our approach the regulariza- 
tion of putative singularities is not achieved with surface 
tension, but by having a minimal size bump, similarly 
to the regularization of DLA. Our rules of growth with 
(3 = 2 and Aq chosen once and for all, guarantee that 
every layer of growth has exactly the same area. This in 
the continuous time Laplacian Growth model translates 
to a particular choice of the time step dt. Clearly, one has 
freedom in choosing dt, or of the size Aq in each layer, as 
long as this does not affect the nature of the growth. In 
particular we can have Aq chosen such that the maximal 
physical bump is of constant area. Once Aq is chosen, the 
sharpest feature that can be achieved is a bump of size 
Ao, and the worst possible "singularity" is a line of such 
bumps, exactly as in DLA. Thus the putative cusp sin- 
gularity of Laplacian Growth is avoided in a manner that 
is identical for all the growth models in our 2-parameter 
family. 



The conformal map (uj) admits a Laurent expan- 
sion 



F 



(n) 



$(")(c^) = Fi'^'cj + Fo'") + —^ + 



(12) 



The coefficient of the linear term is the Laplace radius, 
and was shown to scale like 



where S is the area of the cluster, 

n 



(13) 



(14) 



Note that for /3 = this and equation (||) imply that S = 
uXq. Indeed for (3 = this estimate had been carefully 
analyzed and substantiated (up to a factor) in [|ll|. On 

(n) 

the other hand is given analytically by 



F^''^ = n v/otaI) 



(15) 



/c=l 



and therefore can be determined very accurately. 

The conclusion from the calculations presented in ||l|, ||] 
is that for C > the fractal dimension of the growth pat- 
terns depends continuously on the parameters, growing 
monotonically upon decreasing (3 or increasing C. It is 
quite obvious why increasing C should increase the di- 
mension. By forbidding particles to overlap we simply 
force them into the fjords, not allowing them to hit the 
tips only (as is highly probable). Also decreasing /3 in- 
creases the dimension, since we grow larger particles into 
the fjords, whereas increasing (3 reduces the size of parti- 
cles added to fjords and increases the size of particles that 
accrete onto tips. In particular we argued that DLA and 
our discretized Laplacian Growth cannot have the same 
dimensions, putting them in different universality classes. 
In the rest of this paper we make these observations more 
quantitative and precise. 



III. MULTIFRACTAL PROPERTIES 

A. Generalized Dimensions 

The fractal dimension in the (3 — C family of models, 
D{(3, C), is introduced as the exponent relating the area 
of the cluster Sn to its linear scale (which is measured by 
the (dimensionless) Laplace radius F^^^''): 



(16) 



In this equation Aq = Aq^*'^^'^''. The multifractal ex- 
ponents ll^ are defined in analogy to those for DLA in 
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terms of the moments of the (dimensionless) electric field 
E{s) on the bomidary of the cluster |p^ , 

^^(9-1)) _ (^^("))-(«-i)£'<,(/3,C') 

^ (5„/Ao)-(«-i)^'(^'^)/^(^'^), (17) 

where (• • •) represents the harmonic average for the C) 
clusters in question. Note that these exponents are for 
a fixed size partition with boxes of length ■\/Ao, with 
asymptotics for an infinitely large cluster. A supremum 
over arbitrary partitions may lead to different exponents, 
cf. 

This result translates immediately to the multi- 
fractal fluctuations of the bump areas A„ added in the 
mathematical plane. As ~ En^^'''^ , where En is the 
field computed at z{s) = $(")(e'^"). We therefore write 

^ (5„/Ao)-('+««^(^+«'+^(^'^)/^('''^^ (18) 

Specifically we can derive the following important mo- 
ments 



{VK) - (5„/Ao)-(i+'^/^)^'^+'^/^>/^ 

(A„) ~ (5„/Ao)-(2+/3)^3+./i3 (19) 

(Af/(2+/3)) ^ (5„/Ao)-^^^+^/^ 

(20) 

where naturally all the dimensions are functions of C). 
We can also estimate the way in which the maximal bump 
areas scale 

(21) 

Consider now the addition of one layer of p particles 
to the growing cluster. We can rewrite Eq. (]ll|) as 



C= (l/7r)pVA. 



(22) 



where we have introduced the notation An to represent 
the average over a layer of p particles: 



1 ^ 

n ^ — ^ 



(23) 



fc=i 



For our considerations below it is important to relate 
the layer averages A^ to harmonic averages (A^). This 
relationship may very well depend on the value of [3. The 
two case that are of highest interest to us are [3 — Q and 
(3 — 2, and we will examine them separately. 



averages vs. the number of layers for the case q = 1, 
(3 — and four values of C. In Fig. || we show the same 
for the case q — 0.5, P — and the same four values of 
C. Examining the results it appears that for the higher 
values of C we can assume that in the scaling sense 



Xl ~ {XI) , (3 = 



(24) 



Note that for smaller values of C the evidence is not as 
clear cut as for higher values. The number of points p 
in each layer is relatively small and the layer average is 
highly fluctuating. Nevertheless, even for the case C = 
0.01, if we perform a running average on the layer average 
data, we converge very well onto the harmonic average. 
We therefore propose to proceed with the conjecture that 
Eq. (^) is correct for all the values of C and P = 0, and 
investigate the implications of this scaling relation for the 
cases for which it is correct. An immediate consequence 
of Eqs. M) and {M) is that 



P ~ C/( Va;) ^ C(5/Ao)^^/^ 



(25) 



We note that this means that p ^ oo asymptotically for 
every value of C, while p/n — > 0. 
Next observe that by definition 



In light of Eq. ( pj| ) we write 



5, 



n+p 



(n+p 



D 



On the other hand we estimate 



1 + aDpXn 



Sn+p ^ 1 ^ pAp 

and comparing with (^^ we find 
S ~ 

aDXn 

If Eq.(p3) is used, we find finally 



Sn ~ Xq 



2D3/D 



(26) 



(27) 



(28) 



(29) 



(30) 



from which we derive the well known "electrostatic rela- 
tion" 



IV. SCALING RELATIONS FOR THE 
FRACTAL DIMENSION D 

A. The case (3 = and C > 

We examine the relationship between layer and har- 
monic averages numerically. In Fig. we show the two 



D = 2D3. (31) 

This result was known for C = jist , and is generalized 
here, under the conjecture ( p4| ) to all values of C. 

Let us consider now the probability to hit at the point 
of maximal radius. We propose that for any finite C the 
probability for this event is finite. We stress that this 
"point" is actually a region on the interface of size \/Ao 
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FIG. 1: Layer and harmonic averages of A„ as a function of 
the number of layers, for /3 = 0. Panels a-d: C = 0.01, 0.1, 0.3 
and 0.5 respectively. 
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FIG. 2: Layer and harmonic averages of A^^'^ as a function of 
the rmmbor of layers, for /? = 0. Panels a-d: C = 0.01, 0.1, 0.3 
and 0.5 respectively. 
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in every layer. In particular, we expect that the growth 
process will hit the point of maximal radius every finite 
number of layers, where this number is of the order of 
1 jC . We also know for sure that we have at most one hit 
per layer since particles cannot overlap in the dynamics. 
Consider now the scaling of the size of the growth pat- 

fn) 

tern which is measured by . First we know that 
F^^^ - (S'/Ao)i/-°, and therefore 

dF'^'^ldS^iSIX^fl^-^IX,. (32) 

On the other hand, we estimate the same object using the 
following argument: the maximal radius i?^") increases 
by \f\) every time that it is hit. This occurs every 1/C 
layers in which p particles were added. Therefore 

^ ^ (33) 

Comparing Eqs. (^) and (^, using Eq. ( |25| ) we obtain 
the scaling relation 

D = l+D2 (34) 

Using the inequalities between the generalized dimen- 
sions and Eq. ( |3l| ) we write 

D -1 = D2>D3^ D/2 for aU C > C„ , (35) 

which is equivalent to 

D = 2 for aU C > C^r (36) 

In other words, we conclude that along the line (3 = 
in the phase diagram /3 — C, there exists a transition to 
growth patterns of dimension 2. 

Since our arguments are not rigorous and the result 
quite surprising, we will examine the assumptions using 
an additional consideration. From Eqs. ( p4[ ) and ( |35| ) 
follows that I?2 = 1, and from Eq. ( p5| ) it then follows 
that p scales like 

p - S^/^ , C>C„ ■ (37) 

This prediction is examined directly in Fig. ||. We see 
that it is obeyed extremely well for all the values of C > 
0.1, and it is not in contradiction with the data even for 
C — 0.01. We therefore cannot exclude the possibility 
that Ccr = 0. 

To gain intuition to the meaning of this result we show 
in Fig. ^ the actual growth patterns for /3 = and 
C = 0.01, 0.1, 0.3 and 0.5 respectively. To plot these 
figures we find all the exposed branch cuts on the unit 
circles which are associated with the bumps added in the 
growth process (see for details). Then we plot the 
image of all these points under the conformal map and 
connect them by lines. Thus we are guaranteed that what 
is plotted is the actual contour of the growth pattern, of 
the image of the unit circle in the mathematical domain. 




n 



FIG. 3: The number p of bumps in a layer vs. the number 
n of bumps in the growth pattern, in a log-log plot. From 
top to bottom that is shown for C = 0.5, 0.3, 0.1 and 0.01 
respectively. In white lines we show the scaling laws p ~ n^^^; 
this law fits the data for C > 0.1 and is not in contradiction 
with the (noisy) data even for C — 0.01. 



with all the fjords fully resolved. We see that even with 
the lowest value of C the branches appear to gain sub- 
stance as they grow, having a width which is larger than 
VAo (the typical corrugation of the interface). Conse- 
quently it is not impossible that D — 2 even for the 
lowest values of C > 0. If this is so, it is NOT due to the 
existence of an ultraviolet cutoff, but due to the finite- 
ness of C. With C = (the DLA hmit) the serial algo- 
rithm favors strongly truly fractal patterns. The parallel 
growth algorithm with finite C squeezes more substance 
into the fjords, reducing that tendency. For higher val- 
ues of C it becomes obvious that the growth patterns are 
2-dimensional, and for C — 0.5 the pattern grows like a 
roughened disk. The main conclusion of this analysis is 
that we certainly cross somewhere along the line P — 
into growth patterns that are 2-dimensional. Whether 
or not the critical value of C is finite or zero cannot be 
determined by numerics alone. 

If we accept the possibility that even the lowest val- 
ues of C are associated with growth patterns that are 2- 
dimensional, then we should stress that standard ways of 
estimating the dimension of these clusters, especially for 
the lowest value of C, may fail to discover this fact. For 
example, we can compute F^^^ and then, using Eq.(p^, 
attempt to extract the dimension from log-log plots or 
F^""^ against S. This method works very well for the 
fractal case, but it does not appear to do so well for the 
cases at hand. In Fig. || we show such log-log plots for all 
the clusters of Fig. m We see that even with 100 000 par- 
ticles the dimension estimate is way below the suspected 
D — 2, except for C = 0.5. In fact, any practitioner in 
the fractal field would be happy to interpret the scaling 
obtained for C = 0.01 as an indication that it is in the 
same universality class of DLA with dimension very close 
to D = 1.71. While we cannot state confidently that for 
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FIG. 4: Clusters for /3 = 0. Panels a-d: C = 0.01, 0.1, 0.3 and 
0.5 respectively. Note the areas significantly larger than the 
UV-cutoff Ao which appear already for C = 0.01. 




FIG. 5: The first Laurent coefficient F^'"' as a function of the 



area for /? = and C = 0.01, 0.1, 0.3 and 0.5. 
mension D is obtained for the slope via F^"' 



The fractal di- 



C = 0.01 the growth pattern is 2-dimensional, we stress 
that the dimension estimates obtained from log-log plots 
can be only taken as lower bounds on the true dimension, 
and these may not be very sharp. 

A possibly better way to measure the dimension would 
be through the result (Q) when it holds. We have very 
good methods to determine the correlation dimension 
D2, going back to the Grassberger-Procaccia algorithm 
]l6t . To this aim we choose randomly m = 100000 points 
{^i}"=ii and compute their positions on the interface of 
the cluster Zi = $("^(e'^'). Next we compute the corre- 
lation integral 



C(2)(r) = ^e(|z,-z,|-r) , 



(38) 



where Q{x) is the step function, being 1 for x < and 
for a; > 0. The correlation integral is known to scale 
according to 



(39) 



In Fig. g we display this object in a log-log plot as a 
function of r. All the values of C agree with a correlation 
dimension of D2 = 1, as can be seen from the plots at 
small scales. For those values of C for which Eq. (^) is 
correct this leads to the aforementioned result D = 2. 



B. The case /3 = 2 and C > 

The next interesting family of growth patterns that we 
focus on is obtained for (3 = 2 and C > 0, with Laplacian 
Growth expected to be realized for C = 1. We find that 
for /3 > the numerics does not support the scaling re- 
lation (24). In Figs. (0) and (H) we show the layer and 
harmonic averages for /3 = 2, and it is obvious that in 
this case 



(40) 




10000 



radius 



FIG. 6: The correlation dimension D2 for (i — and C — 
0.01,0.1,0.3 and 0.5. The thick hne has slope 1, indicating 
that D2 = 1 and therefore D — 2 for all shown C. 



in the scaling sense. 



Once we have lost the scaling relation (24) we cannot 
argue that D — 2 ior any value of C > 0. We will find 
numerically that along the line (3 = 2 we indeed find frac- 
tal patterns, (and cf. the next section); nevertheless even 
along this line there exists a transition to 2-dimensional 
patterns, albeit at a finite and rather high value of C. 
Next we want to estimate this value. 



V. CONJECTURE: LAPLACIAN GROWTH IS 
2-DIMENSIONAL 

In this section we motivate our conjecture that Lapla- 
cian Growth patterns are not fractal patterns at all, but 
rather patterns of dimension 2. We have to be a bit cir- 
cumvent, since as explained in we cannot directly run 
our algorithm for the (3 — C model for values of C higher 
than about 0.65. The reason is that it becomes impossi- 
ble to fill up, by random selection of points on the unit 
circle, a full layer of bumps on the physical interface. 
Therefore our aim is to find a line in the [3 — C phase di- 
agram that separates fractal D < 2 from 2-dimensional 
growth patterns. That such a line must exist we can con- 
vince ourselves by examining the family of growth models 
that are seen for /3 = —1, see Fig. ^ Obviously these are 
2-dimensional. The family of growth patterns obtained 
for /3 = were shown in Fig. ^, and as we said above, 
there must be a cross over 2-dimensional patterns in this 
family. Going up to /3 = 1 we show the growth patterns 
in Fig. |ll|. In this case the images indicate that for the 
lower values of C the growth patterns are fractal, whereas 
for higher values of C they become 2-dimensional. Thus 
the line of separation that we seek in the P — C phase 
diagram appears to cut the P = 1 line. Finally, in Fig. 
|l2| we present the family of growth patterns obtained for 
(3 = 2. It appears that the transition line intersects also 
the /3 = 2 line. 

AU the patterns exhibited in Figs. |, are grown 
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and 0.5 respectively. 
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and 0.5 respectively. 




FIG. 9: The first Laurent coefficient _F^'"' as a function of the 
area for /3 = 2 and C = 0.01, 0.1, 0.3 and 0.5. The fractal di- 
mension D is obtained for the slope via 



with a fixed size Aq. Consequently, for /? > the actual 
mean size of the bumps in the physical space decreases as 
the cluster grows, while it increases for /3 < 0. This may 
lead to worries, i.e. that for /3 > the growth arrests and 
that for /3 < the increase in the size of the bumps leads 
to coverage of fjords, such that the 2-dimensional pat- 
terns shown in Fig. would be an artifact. To disperse 
these worries we have considered alternative growth al- 
gorithms with varying the size of Aq . The first such algo- 
rithm is obtained by requiring that the total area covered 
in each layer of growth is constant, i.e. 



j]A„+fc|$'(")(e^^"+'=)|2 
fc=i 

p 

= Ao(n)^|$'("ne*''"+')r'^ = C . (41) 
fc=i 

Note that for constant coverage C this rule coincides with 
fixed values of Aq for /3 = 2 (cf. Eq.(pT|)). In the second 
algorithm we choose the maximal size of the bump in the 
physical plane to be constant from layer to layer: 

Ao(n) max{|$'(")(e*^"+'=)r''} = C . (42) 

k—l 

This rule coincides with fixed values of Aq for /3 = 0. 
We found that in all cases the patterns shown above re- 
main invariant to the change of the algorithms. Thus we 
submit that the figures shown can be fully trusted. 

To find the line that separates fractal from 2- 
dimensional patterns we estimate the dimensions directly 
from log-log plots of F*^") vs. S. We have seen above that 
such estimates are lower bounds to the actual asymptotic 
dimension. As these logarithmic plots are invariably con- 
cave, we can use the slope at the largest values of area 
available as a measure for the lower bound on the dimen- 
sion. In Fig. O we show the three lines obtained by 
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FIG. 11: Growth patterns for /3 = 1. Panels a-d: C = 
0.01,0.1,0.3 and 0.5 respectively. 
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FIG. 12: Growth patterns for /3 = 2. Panels a^d: C 
0.01,0.1,0.3 and 0.5 respectively. 
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FIG. 13: The phase diagram in the — plane. The data 
points in crosses, triangles and circles represent values of (3 
and C for which the radius-area relationship predicts D = 
1.90, D = 1.95 and D = 1.99 respectively. The lines are 
quadratic fits. We propose that the region below the lines 
represents 2-dimensional growth patterns, and see text for 
details. 



searching, for a given value of C, the value of (3 for which 
for the first time the dimension estimated from i^^") vs. 
S crosses the value D = 1.90 (upper curve), D = 1.95 
(middle line) and D = 1.99 (lower curve). We propose 
that the last two Unes may very well be already beyond 
the true line that separates fractal from D — 2 asymp- 
totic dimension. From the discussion of Sect. IV A we 
cannot even exclude the possibility that the transition 
line obfuscates the /3 = line. All the region below the 
lower line is almost surely representing patterns of D = 2, 
but we strongly believe that this is the case also for the 
middle line. The lines were obtained by finding, as ex- 
plained, the values of f3 yielding D=1.90, 1.95 and 1.99 
respectively, and then fitting to the points a quadratic 
function. Next we extrapolated the three fits to values 
of C that arc not readily available in our algorithm. The 
three fit lines intersect the ft = 2 line at C = 0.73, 0.78 
and C = 0.79 respectively. Wc thus propose that the 
value C = 1 for /3 = 2 is comfortably within the region of 
2-dimensional patterns in this phase diagram. 



VI. CONCLUSIONS 

We have presented a careful numerical study of a 2- 
parameter model of growth patterns that generalizes and 
interpolates between Diffusion Limited Aggregation and 
Laplacian Growth Patterns. The model gives rise to a 
rich plethora of growth patterns, with fractal dimensions 
that depend on the values of the parameters jS and C. 
For /3 = and C = we obtain DLA. Laplacian Growth 
patterns have (3 = 2 and C = 1, but we cannot probe 
the value C = 1 within our algorithm. Since our aim, 
in part, is to demonstrate that Laplacian Growth pat- 
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terns are not fractal, we resorted to examining the phase 
diagram (3 — C. We estabUshed, on the basis of scal- 
ing arguments, simulations and visual observations, that 
this phase diagram contains a line of transition between 
fractal and 2-dimensional growth patterns. Wc have es- 
timated the position of this line, and demonstrated that 
Laplacian Growth patterns belong safely in the region of 
2-dimensional growth patterns. 

One should point out that the statement that Lapla- 
cian Growth are 2-dimensional does not mean that it is 
a growing disk. To the eye the patterns can look frac- 
tal, and in fact radius-area log-log plots might initially 
even indicate that the dimension is low, and maybe of 
the order of the dimension of DLA. Deep fjords may ex- 
ist in the structure. The relevant question is whether 
the growing branches of the structure contain substance 
(area) and whether this area is growing relatively with 
the growth of the pattern. The growth pattern shown 
in Fig. 11 panel d is a case in point. It looks fractal to 
the naked eye, but careful examination shows that the 
branches have area. Thus one needs to decide whether 
this area is due to some ultraviolet cutoff length, or does 
it grow systematically beyond what is expected on the 
basis of the existence of such a cutoff. 

Before closing we reiterate that our demonstration that 



Laplacian Growth patterns are 2-dimensional is not di- 
rect. We cannot, within our algorithm, grow C = 1 pat- 
terns. We therefore leave this at the moment as a con- 
jecture. It remains a theoretical challenge to show that 
this conjecture is indeed provable by direct mathemati- 
cal analysis. We also leave for future work the question 
whether the /3 = line represents 2-dimensional growth 
patterns for all C > 0. Finally we propose that future 
work may make use of the fractal patterns along the line 
C = 0, /3 > for further fundamental studies of DLA and 
related phenomena. 
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